Mapping from Architecture to Dynamics: A Unified View of Dynamical 

Processes on Networks 
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Although it is unambiguously agreed that structure plays a fundamental role in shaping the dy- 
namics of complex systems, this intricate relationship still remains unclear. We investigate a general 
computational transformation by which we can map the network topology directly to the dynamical 
patterns emergent on it — independent of the nature of the dynamical process. We find that many 
seemingly diverse dynamical processes such as coupled oscillators and diffusion phenomena can all 
be understood and unified through this same procedure. Using the multiscale complexity measure 
derived form the structure-dynamics transformation, we find that the topological features like hier- 
archy, heterogeneity and modularity all result in higher complexity. This result suggests a universal 
principle: it is the desire for functional diversity that drives the evolution of network architecture. 
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Recent advances in the realm of complex networks have 
furnished us with a new paradigm to understand and 
characterize complex systems [H, H), H| ; from the human 
brain that is composed of billions of interconnected neu- 
rons, to our society with six billion cooperating individu- 
als. The discovery of the distinctive scale free [J] and the 
small world [5| structures has fundamentally altered our 
view of these networks. A variety of other complex topo- 
logical features, such as high clustering [||, hierarchical 
ordering Q, and degree mixing 0, are also emerging as 
important to the overall behaviour of network systems. 
However, recent progress has focussed mainly on the un- 
derlying topology [lj , the effort to unravel the system's 
dynamics or function has been less advanced 0, 0, HJ- 
With the increasing capability to capture simultaneously 
the time dependent activity of many components in com- 
plex systems OH (such as the multiple electrode record- 
ings and gene expression patterns), unraveling the intri- 
cate relationship between the structural and functional 
characteristics has become a problem of utmost impor- 
tance, and hints toward generic organizing principles and 
a deeper understanding of "complexity" . Generally, the 
topological descriptors fail to explicitly capture the dy- 
namical aspects. In order to characterize the dynamics, 
one must implement the particular dynamical process on 
networks via extensive numerical simulations. 

In this paper we propose a methodology that allows 
one to predict the collective dynamics (or functions) of 
a complex system directly from the underlying topol- 
ogy. Yet, this methodology is independent of the details 
of dynamical process. This is achieved by constructing 
node interaction profiles through a kernel function, which 
quantitatively identify the role of each node in shaping 
the integrated dynamics and thus captures the way nodes 
are dynamically interacting with one another. Notably, 



we find that synchronization of coupled oscillator, en- 
semble neuron firing, epidemic spreading and diffusion 
process can all be unified under this theoretical picture, 
suggesting that some universal mechanisms may govern 
the overall dynamical behavior of the seemingly diverse 
complex systems. Based on this transformation we pro- 
pose a "function-driven" multiscale complexity measure 
by virtue of the adjustable kernel bandwidth, which un- 
ravels the functional organization of a network at differ- 
ent levels. We find that structural complexity measures 
such as topological heterogeneity [l[ , modularity @ , hier- 
archy Q, and nontrivial correlation Q all translate into 
higher functional complexity, indicating that the need for 
multiple function governs the structural evolution of net- 
works. 

Consider a complex system of N coupled dynami- 
cal units, whose equations are described by fill ]: x\ = 
F(xi) + (T^ = i AijH(Xj),i = 1,2, ...,7V where x l = 
F(xi) governs the dynamics of each component, H is a 
fixed output function, a represents the coupling strength, 
and A is the coupling or adjacency matrix. The prob- 
lem is now, given the topology A of the network, can we 
infer qualitatively the dynamics without implementing 
and computationally simulating it? If so, what does this 
tell us about the behavior of different dynamical units 
within the same network structure? By "dynamics" we 
refer to the pattern of temporal correlations among out- 
puts #i(t), i = 1, N. This has been intensively studied 
in chaos community and in brain researc h g roups, un- 
der either the banner o f sy nchronization [ill Hil [l3| or 
functional connectivity 14j. 

Obviously, to predict the collective dynamics, we must 
identify the effective influence of each unit to others, and 
distinguish their specific, functional roles in shaping the 
dynamics of each unit. A fundamental feature of com- 
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plex interacting systems is the presence of interactions 
across all scales. To describe the fact that directly cou- 
pled components usually exert a stronger influence on 
each other and the impact attenuation through inter- 
mediaries nodes, we adopt a simple, monotonically de- 
creasing function known as a "kernel" K. For a pair 
of nodes i and j, we then define their effective inter- 
action Rij , in the form of K (xij , h) , where X % j IS their 
shortest distance, K is a non-negative, symmetric kernel 
function satisfying J R K(x)d(x) = 1, J R xK(x)d(x) = 0, 
lim x _ >oo K(x)d(x) — 0, and h is the bandwidth that con- 
trols the width of the kernel. For Gaussian kernel, the 
interaction matrix R will read: = exp(— x^/2/i 2 ). 

Note that the i-th row of -R, Ri portrays the effective 
interaction node i receives from all its neighbors, which 
we call interaction profile of node i. This vectorial pro- 
file systematically encodes the identity of all neighbors of 
node i as distinct driving forces (with different intensities 
determined by the kernel) to its own dynamics, and does 
not depend on specific choice of kernels. Therefore Ri 
defines the unique "status" of node i. To further predict 
the dynamical correlation or functional connectivity Fij 
between node i and j, we can calculate the similarity be- 
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tween their interaction profiles Ri and Rj : = rffiTRn 
The similarity Fij provides a unique clue to evaluate 
the dynamical proximity between the components. Unit 
i and j subject to a large number of common inputs (up 
to higher orders) are more likely to behave similarly. In 
this case, their profiles Ri and Rj will largely coincide by 
sharing many common entries, leading to a large i*y - 
approaching 1. Conversely, a pair of units with few com- 
mon drives tend to be independent and thus have a Fij 
near 0. A great advantage of this "kernel" formalism lies 
in the adjustable bandwidth h which can evaluate differ- 
ent levels of function of the network at various topological 
scales. In the following we will demonstrate how the col- 
lective dynamics of various dynamical processes can be 
predicted using the above mapping from structure A to 
dynamics (or function) F. 

We start with synchronization phenomena, which are 
widely observed in nature and occupy a privileged posi- 
tion in understanding collective behavior in various dis- 
ciplines [l2| . Recently the interplay between a network's 
structure and its synchronization dynamics has attracted 
significant attention [3, EH- Here we use Kuramoto 
model defined on various networks as a prototype exam- 
ple. It is govern by §i = 0Ji + Y!rjLi sm(8j -#;),£ = 
1,2, ...,7V, where LOi is the frequency of phase oscilla- 
tors (uniformly distributed in (0,1)), a is the coupling 
strength, A is the adjacency matrix, and k is the mean 
degree. Specifically, we will show that obtained at 
different kernel bandwidth h can provide a good predic- 
tion of how the collective dynamics evolves with a. 

With a small coupling a, the oscillators are mostly 
independent. The dynamical distance between outputs 
of node i and j, defined as Z)y = (6i(t) — 9j(t)) (8(t) 




FIG. 1: The correlation coefficient p between D t j and 1 — Fy 
for coupled phase oscillator on (a) random (ER) network with 
500 nodes (mean degree is 5) and (b) a collaboration network 
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FIG. 2: Visualization of matrixes (a) Dij and (b) 1 — Fij ob- 
tained form Fig. Q](b), i-e., phase oscillators on collaboration 
network, with a = 0.6 in computing D„ and h — 2.6 for Fij. 

are wrapped to [0, 2-7r] and (■) means time average) will 
be non-zero and constitute a narrow distribution. When 
a is large, the whole network reaches complete synchro- 
nization, and distribution will be a narrow peak again 
near 0. For intermediate a, various functional clusters are 
formed, with distribution broadening. We find that 
the distributions at various a are exactly reproduced 
by F^ using different h. At a small h, all entries in Ri 
are almost except the ith, meaning each node only has 
impact on itself. The RiS are mostly orthogonal, thus Fij 
will centralize at 0. By contrast, the kernel becomes flat 
at a large h, leading every node to exert similar influence 
on all others. The Fij then concentrates near 1 as all RiS 
are almost identical. For medium h, the distribution of 
F^ broadens within [0, 1] with the peaks corresponding 
to the formed functional clusters. 

These examples show that the kernel bandwidth h is 
playing a role directly analogous to coupling strength a, 
and F^ offers a good prediction of the collective dynam- 
ics D^. To further verify this, we first get D^ by imple- 
ment phase oscillators on networks (a random network 
and a modular collaboration network [l6| are used here 
as examples). We then plot the correlation coefficient p 
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between D^ and 1 — Fij that is obtained at various h, see 
Fig. Q] (we use 1 — F^ because it measures dissimilarity 
similar to Dy, while F^ is a similarity measure). Wc 
find that for a given a, there is always an optimal kernel 
bandwidth h max that attains a maximum similarity be- 
tween Dij and 1 



F^ , with h max being proportional to 



a. The matrices and 1 — Fij demonstrate vary simi- 
lar patterns (see Fig. [2|), reflected by a large correlation 
coefficient p. We find that p generally takes a large value 
for medium and high coupling a, where oscillators have 
self-organized into functional clusters. For weak er, the 
oscillators are largely independent, thus is somewhat 
random and cannot be accurately fitted by F^ . 

Now we turn to a concrete example in neuroscience, 
the coupled neural oscillators in cortex, which commu- 
nicate by non-smooth, pulse- like firings. The population 
dynamics of the neurons, like the synchronous firing plays 
a vital role in cognitive function of the brain. Therefore 
understanding how connectivity patterns influence the 
emergent dynamics is of special concern in neuroscience. 
Here we try to approach the population neuronal dy- 
namics directly from the underlying anatomy. In partic- 
ular, we couple the FitzHugh- Nagumo neurons through 
real networks (by excitatory synapse with synaptic con- 
ductance g) to get and check if it can be predicted 
by Fij. We use two networks possessing key topologi- 
cal properties of the cortex, i.e., small- world, hierarchi- 
cal and modular structure. One is the neural network 
of Caenorhabditis elegans whose anatomy has been iden- 
tified. The other is a hierarchically organized modular 
network Il3| with two hierarchical levels. 




FIG. 3: The correlation coefficient between Dij and 1 — Fij 
for (a) C. elegans network with 297 neurons and each has 14 
synaptic couplings on average, (b) Hierarchical network with 
480 nodes [131 ] - Each node has 20 links to the most internal 
community (formed by 30 nodes) , 2 links to the most external 
community (120 nodes that form four 30-nodes groups), and 
1 more link to any other node. See supplement material for 
the detailed parameters of the FHN neurons. 

The synaptic conductance g determines the amplitude 
of pulse conducted to post-synaptic neurons. It plays the 
same role as a in coupled phase oscillators. The neurons 
fire almost randomly with a small g, and begin to form 



synchronous firing as g increases, giving rise to coherent 
oscillations. We define the collective dynamics of neuron 
firing as Aj = (\\fi(n) - fj(n)\\), where fr(n) is the num- 
ber of firings within time window n for neuron i. We find 
that the functional similarity F^ obtained purely from 
the network structure shows a high correlation with Dij 
for both the C. elegans and the hierarchical network, see 
Fig. Therefore we can precisely predict the dynamics 
of neuronal populations. Moreover, the kernel bandwidth 
h plays a role similar to synaptic conductance g. For a 
given g, there is always an optimal bandwidth h max un- 
der which F^ best fits D^ , and this h max is proportional 
to g. 

The above mapping can predict the collective behavior 
not only for the coupled dynamical systems, but also for 
the general diffusion process on networks like epidemic 
spreading |l8|. The collective behavior here means the 
correlation among the epidemic dynamics of the indi- 
viduals, with the dynamics of node i being a discrete- 
time stochastic process Si(t) (1 indicates infected and 
for healthy). The collective dynamics is then defined as 
D^ = (\\S'i(n) — S'j(n')\\), where 5-(n) is coarse grained 
from Si(t) by counting the number of 1 in time window 
n. We find that is again nicely predicted by Fij un- 
der a suitable h, see Fig. 21 where we run SIS models on 
two typical networks. Interestingly, we find h max relates 
closely to the effective spreading rate A = v/6, where v 
and S are infection and recovery rate of an individual. In 
Fig. [4] we see that under small and large A is best 
fitted by a kernel with small and large h, respectively. 
This is because an infected node recovers quickly at a 
small A and has a small influence range. A large A makes 
the node persistently infective to even its higher order 
neighbors, thus is described better by a wider kernel. 
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FIG. 4: Correlation coefficient between Dij and 1 — Fij for 
epidemic spreading on (a) scale free BA network (500 nodes 
with mean degree 5) and (b) email network jl9( (1133 nodes). 



Having established the mapping from structure to 
function, we are naturally led to the fundamental prob- 
lems: why do the various structures such as modular- 
ity, hierarchy and degree mixing exist in real networks? 
What are their roles in shaping the function? Here we 
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provide an explanation under the framework of complex- 
ity by exploiting the above structure-function transfor- 
mation. Understanding complexity 2(| 21] has long been 
a grand challenge that spans a wide variety of fields. 
Sporns et. al. proposed "neural complexity" to mea- 
sure functional connectivity by implementing Gaussian 
dynamics 22| on networks. Evolutionary speaking, the 
survival of a complex system hinges crucially on the ver- 
satility of its function. Inspired by this, we propose to 
evaluate network complexity by the dynamical patterns it 
can support. Recall that in coupled oscillators, the evolu- 
tion of the collective dynamics versus coupling strength 
are well captured by under various h, thus we can 
approach the function directly through F^ . Here we use 
entropy, defined as H{F) = Y^T ~ Pi l°g(Pi)> to charac- 
terize the diversity of Fij obtained at different h. We 
call it multiscale entropy, as it portrays the functional 
complexity of the network at various topological scales. 
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FIG. 5: Multiscale entropy for (a) scale free BA, random ER 
network, assortively and disassortatively mixed BA network. 
All networks have 1000 nodes and mean degree 20. (b) Hi- 
erarchical networks (used in Fig. [3} and (c) Collaboration 
network [l6| ]. Here we normalize H by a factor H m — log(m), 
which is the entropy for uniform distribution, and m is the 
bin number. 

Now we examine how the intricate topologies can in- 
fluence the network complexity by computing multiscale 
entropy. As can be seen in Fig. Ufa), the degree hetero- 
geneity (scale free distribution) and degree mixing Q (as- 
sortive and disassortative) both lead to higher H at var- 
ious scales, as they promote differentiation or facilitate 
the formation of modules that can engage into various 
functions. Remarkably, we find that hierarchical modu- 
lar networks, which are widely observed in biological and 
social systems 0, demonstrate notably higher complex- 
ity than their random counterparts, see Fig. [SJb) and (c). 
This is because the multiscale modular structure can pro- 
vide different levels of function persistently at different 
scales. The fact that all these distinctive topologies lead 
to higher functional complexity provides important in- 
sights into the evolutionary mechanism of real networks. 
This suggests a common principle: the demand for func- 
tional capability shapes the network architecture during 
the development of a physical network. 



In summary we have introduced a novel approach that 
can map the topological structure of a network directly 
to its functional organization independent of the details 
of dynamical processes. Our method not only provides a 
good prediction of the collective network dynamics, but 
also allows us to conceptualize the "complexity" of net- 
works through the emergent functions conveniently. This 
mapping can furthermore be considered as a promising 
scheme to other problems like inverse engineering and 
controllability. The explicit relation between structure 
and function will provide unique clues to infer the struc- 
ture back from dynamical patterns, possibly with con- 
straints like sparseness of connectivity. The control over 
the general networked systems can also be expected to en- 
hance conveniently by locating the most sensitive nodes 
or links, the removal or rewiring of which would result in 
better network performance. 
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